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Abstract 

D . 

An extended version of 4-d SU(2) lattice gauge theory is considered 
in which different inverse coupling parameters are used, Ph = '^I 9h 
plaquettes which are purely spacelike, and Pv for those which involve 
the Euclidean timelike direction. It is shown that when I3h = oo the 
partition function becomes, in the Coulomb Gauge, exactly that of a 
set of non-interacting 3-d 0(4) classical Heisenberg models. Long range 
order at low temperatures (weak coupling) has been rigorously proven 
for this model. It is shown that the correlation function demonstrating 
spontaneous magnetization in the ferromagnetic phase is a continuous 
function of gn at gn = and therefore that the spontaneously broken 
phase enters the {Ph, Pv) phase plane (no step discontinuity at the 
edge). Once the phase transition line has entered, it can only exit at 
another identified edge, which requires the SU(2) gauge theory within 
also to have a phase transition at finite /?. A phase exhibiting spon- 
taneous breaking of the remnant symmetry left after Coulomb gauge 
fixing, the relevant symmetry here, is non-confining. Easy extension to 
the SU(N) case implies that the continuum limit of zero-temperature 
4-d SU(N) lattice gauge theories is not confining, in other words gluons 
by themselves do not produce a confinement. 
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order 



For many years it has been believed that SU(N) lattice gauge theories 
exist in a single phase. Since confinement is a property of all compact lat- 
tice gauge theories at strong coupling, this assumption has the consequence 
that these theories also have a confining force (linear potential) in the weak 
coupling continuum limit. This fits rather nicely with phenomenological ev- 
idence that the strong interactions confine quarks. The numerical evidence 
for lack of a phase transition in SU(2) comes mostly from the apparently 
smooth behavior of the specific heat as a function of the coupling parameter, 
(3. However, if the critical exponent v is large enough to produce a large 
negative specific heat exponent a {a = 2 — dv), then the first infinite singu- 
larity may be in a high derivative, and not easily visible numerically. This is, 
for instance, often the case for 3-d spin glasses which are close to their lower 
critical dimension. The argument in favor of the existence a phase transition 
rests on the similarity of lattice gauge theories to ferromagnetic spin models, 
all of which in three or more dimensions, at least for short range interactions, 
have two phases, with the weak coupling phase exhibiting spontaneous sym- 
metry breaking and long range order. Abelian and non-abelian spin models 
differ in detail, but all have phase transitions. A major difference between 
gauge models and spin models are that the former have local symmetries and 
the latter global. This difference can sometimes be erased through gauge 
fixing, however. For instance the 2-d lattice gauge theory is exactly equal 
to a set of non-interacting 1-d classical spin chains in the axial gauge, ne- 
glecting possible boundary effects. Setting all of the links in one direction 
equal to unity using the gauge freedom causes the 4-link plaquette action to 
collapse to a two-link dot product between the remaining gauge links which 
can be reinterpreted as spins \T\ . Another exact mapping is between the 3-d 
Ising gauge theory in axial gauge with the 3-d Ising model For 4-d theo- 
ries, the Coulomb gauge appears to be the most "spin-like." Here the gauge 
freedom is used to maximize the traces of all of the links lying along three 
of the four lattice directions. In this gauge the fourth direction lying links 
become observables and act much like spins, in that they transform under 
a remnant 3-d global symmetry on each spacelike hyperlayer that remains 
after Coulomb gauge fixing. This symmetry can break spontaneously if the 
average of fourth-direction links over a hyperlayer acquires a magnetization. 
It has been shown that if this symmetry does break spontaneously, then the 
corresponding ferromagnetic phase is non-confining, so the average of these 
4-th direction links serves as a local order parameter for confinement [3] . Ac- 
cording to the standard hypothesis, the theory remains in the paramagnetic 
phase for all couplings in the non-abelian case, and undergoes a magnetic 
phase transition only in the abelian case, such as U(l). However, below it 
will be shown that this inconsistent with the known long-range order for 
classical 3-d Heisenberg models, to which the 4-d lattice gauge theories are 
closely related. 

Numerical evidence for a ferromagnetic phase transition based on the 
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behavior of the Coulomb gauge spin-hke order parameter described above 
has been given for SU(2) with the Wilson action. Finite size scaling shows 
an infinite-lattice transition around /3 = 3.2 with a correlation length criti- 
cal exponent v = 1.7 it 0.2[4J. Other simulations [1| which supplement the 
Wilson action with an infinite chemical potential for gauge-invariant S0(3)- 
Z2 monopoles [5], together with a positive plaquette constraint, appear to 
remain in the ferromagnetic phase for all couplings. This would seem to 
indicate that the normal confinement seen with the Wilson action is due to 
strong-coupling lattice artifacts, not unlike the U(l) case. In the following, 
the analytic case for the existence of this phase transition is explored. An 
analytic proof for the opposite hypothesis, namely the lack of a phase tran- 
sition with confinement persisting to the continuum limit, was presented by 
Tomboulis some time ago[6] and updated more recently [7]. This proof is 
fairly complex and some possible flaws have been noted[8]. When an hy- 
pothesis is difficult to prove, it is sometimes worthwhile to attempt to prove 
the opposite, which is the approach taken below. 

The 4-d SU(2) theory can be extended into a larger coupling space by al- 
lowing the coupling for purely spacelike (horizontal) plaquettes to differ from 
that of plaquettes which include the Euclidean timelike direction (vertical). 
The action is 

(2 3 
/5h E E (1 - ^tr(U5,iU5+.-,jUt^^-^,ut p) (1) 
i=\ J=i+1 

E(l - ^tr(Us,iUa+j,4U^ .4,4))] 

Where the C/^ j are SU(2) valued gauge links based at lattice site n in di- 
rection J. The normal Wilson action is simply the fiy = 13 h case, however 
because the renormalization group gives a relation between (3 and the phys- 
ical lattice spacing through the running coupling, the j3y 7^ (in theory can 
also be considered an SU(2) LGT with unequal lattice spacings in the space- 
like and Euclidean timelike directions. So at least in the phase connected to 
the continuum limit (neighborhood of {(in-, /3y) — > (00,00)) the entire phase 
plane excluding the boundaries can be considered an SU(2) LGT. 

Consider now the /3// = 00 theory. This can be seen to be equivalent 
to a set of non- interacting 3-d 0(4) classical Heisenberg models as follows. 
At (3h = 00 all of the horizontal plaquettes will be forced to their largest 
possible value of unity. One can then find a gauge in which all horizontal 
links are unity also, as follows. Set a maximal tree of links to unity on 
each spacelike hypersurface. This is a partial axial gauge in which the final 
trunk of the tree along the Euclidean timelike direction is not completed. 
For instance for an lattice, on a given spacelike hyperlayer all 3-direction 
links can be set to unity except when z = 0. On the 2 = plane all 2- 
direction links are set to unity except for when y = and along the line 
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[y = 0, z = 0) all 1-dircction links except for when x = can be set to unity 
by gauge transformations. Looking at the z = plane, there is a set of x-y 
plaquettes extending backward from a; = (through the periodic boundary 
condition) which have three links set to unity. The plaquette being unity due 
to I^H = oo ensures that the fourth link in the plaquette is also unity. Now 
the same is true for the next row of plaquettes etc., forcing all links to unity 
except for the last row pointing in the positive direction from x=0. These 
are now equal to the gauge invariant Polyakov loop for that direction, and 
all of these links are equal. This same procedure can be extended from the 
{x = 0, y = 0) plane along the z direction to show that the links out of the 
plane are also unity, except for one at the edge along each lattice direction 
which are equal to their neighboring links pointing in that direction. Finally, 
a Polyakov loop symmetry transformation (also a symmetry of the action) 
can be employed to bring these final set of links to unity. In this restricted 
sector, the Polyakov loop symmetry is the whole SU(2) group rather than 
just the center. One could also formulate the theory with open boundary 
conditions, in which case this final step would be unnecessary. Although 
the language of axial gauge was used above, one can see that actually the 
condition for Coulomb gauge has also been met, that the sum of traces of 
all horizontal links is maximized. Coulomb gauge leaves a layered remnant 
symmetry, one global SU(2) per hyperlayer, unfixed. Because Coulomb 
gauge was being sought is why the axial tree above was left uncompleted. 

Since all of the horizontal plaquettes are unity, the vertical plaquettes 
simplify to 

^tr(uJ^4U5+j,4) (2) 
where j, of course, can be 1, 2 or 3. Writing each U as 

3 

U = SQ + i^ SkTk (3) 

k=l 

where the are the Pauli matrices, one can associate an 0(4) unit vector 
s = (sq, Si, S2, S3) with the SU(2) valued link. In terms of the s's it is easily 
verified that the vertical plaquette above is simply the nearest neighbor 0(4) 
dot product, 

-tr(ut 4Ug+j,4) = ss • sg+j (4) 

The action on each hyperlayer is that of the 3-d 0(4) Heisenberg model with 
coupling parameter (inverse temperature) Pv In a gauge theory there are 
never any direct interactions between links longitudinally, and the freezing 
of horizontal links prevents any indirect interactions among hyperlayers, so 
these all become independent Heisenberg models. Each of these has its own 
SU(2) global symmetry, the remnant symmetry from the Coulomb gauge. 
The reason behind this is that an SU(2) gauge transformation which is 
global on the hyperlayer transforms horizontal links in such a way as to 
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leave the trace of each hnk unchanged. Thus the gauge condition, which 
is to maximize these traces, is not disturbed. Therefore, at the level of 
the partition function, the fin = oo SU(2) lattice gauge theory becomes a 
set of non-interacting 3-d 0(4) classical Heisenberg Models. Durhus and 
Frohlich[9^ earlier pointed out a connection between SU(2) lattice gauge 
theory and the 3-d 0(4) model, but did not consider the case of split vertical 
and horizontal couplings which allows for an exact mapping. 

Long range order (LRO) in the classical Heisenberg model has been 
rigorously proven p!0]. This means that a ferromagnetic phase must exist 
at a finite weak coupling (large /3y). Because this is a symmetry broken 
phase, it must be separated from the strong-coupling paramagnetic phase 
by a phase transition, which has been convincingly found by Monte Carlo 
simulation at a coupling /3 = 0.9360(1) [Hj. The LRO proof does not depend 
on the specific symmetry group. Indeed, the existence of phase transitions 
in all ferromagnetic spin models is well-established and non-controversial 
from both analytic and numerical perspectives. The ferromagnetic order for 
the Heisenberg model is also quite robust to the addition other interactions. 
For instance it is still preserved at zero temperature even if up to 20% of the 
interactions are switched to anti-ferromagnetic (beyond this level it enters 
a spin glass phase) [12], and at higher temperatures with smaller but still 
substantial contaminations. Below, it will be argued that the ferromagnetic 
phase, and therefore also the phase transition, continues to exist for non- 
infinite I3h as well. 

In Fig. 1 the fin - phase plane is shown, with the 3-d 0(4) model 
existing on the top edge {Ph = oo). The normal 4-d SU(2) model exists on 
the f3v = Ph line. However, as stated before, the /3y ^ /3h cases can also be 
considered to be SU(2) lattice gauge theories, with a different lattice spacing 
in the fourth direction than the other three. Thus the entire interior, at least 
in the vicinity of the continuum limit in the upper right corner, is an SU(2) 
lattice gauge theory. 

The bottom edge, with fSn = 0, has only the vertical plaquettes left. 
In an axial gauge with the fourth direction links set to unity, the system 
becomes a set of 1-d spin chains, which are paramagnetic except at zero 
temperature {/3v = oo) where there is a phase transition to the ferromagnetic 
ground state. The right edge of the phase diagram, /Sy = oo, is an odd phase 
where all of the vertical plaquettes are forced to unity. Staying with the axial 
gauge, one can see that the 3-d SU(2) theories on spacelike hypersurfaces 
but different Euclidean times are locked to each other, in other words they 
can only fiuctuate in lockstep. This makes for a theory with a 4-d energy, 
but non-extensive 3-d entropy, so it is basically stuck in the classical ground 
state. Both top and bottom ends of this line are ferromagnetic; there is no 
reason for the whole line not to be also. Thus ferromagnetism exists on the 
upper border for f3v > 0.936 and along the entire right-side border. In the 
following, two arguments are given which show that the ferromagnetic phase 
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Figure 1: Possible phase diagram for SU(2) lattice gauge theory with dif- 
ferent horizontal and vertical couplings. Top horizontal axis is 3-d 0(4) 
Heisenberg model and the star shows its known ferromagnetic phase transi- 
tion [IT] . Square and triangle show locations of phase transitions at = 20 
and fSn = found from Monte Carlo simulations using the Coulomb gauge 
magnetization [4|. Diamond is possible spin-glass to paramagnetic transi- 
tion found previously using the real replica method [I3j. Bold regions of 
upper and right edges show portions of the border known to be ferromag- 
netic. Lines drawn are hypothetical phase boundaries guided by the Monte 
Carlo results. 

persists as one enters the interior of the phase diagram. For a symmetry 
breaking phase transition, once the phase-transition line has entered the 
phase plane it must continue to another edge. This is because the order 
parameter, due to the realized symmetry, is exactly zero in the paramagnetic 
region on the infinite lattice, and of course is nonzero in the ferromagnetic 
phase. An analytic function zero in a finite region is zero everywhere, so a 
line of non-analyticity must separate the two phases completely. The only 
place where the phase transition line could terminate is at the only other 
phase transition on the border, namely the lower right corner as shown. 

Fig. 1 also shows the ferromagnetic phase transitions found with Monte 
Carlo simulations using the Coulomb Gauge magnetization as the order 
parameter. One of these was for = 20 where the phase transition on 
the infinite lattice was determined to be /3vc = 1-01(2). This transition 
closely resembles that of the Heisenberg model in that it has similar critical 
exponents and as such is a strong indication that the phase transition does 
enter the interior of the phase diagram. Below, it will be argued that that 
must be the case from an analytic perspective as well. In order for the SU(2) 
lattice gauge theory in the interior of the phase diagram to avoid the phase 
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transition as in the conventional hypothesis, the transition would need to 
be somehow prevented from entering the interior of the phase diagram at 
all. This means that the small terms that arise in the Hamiltonian when 
Ph is backed off from oo would have to destroy the ferromagnetic order, no 
matter how small the coupling {jSu = 4:/g'jj). This seems odd considering 
the known robustness of ferromagnetism in the Heisenberg model. Say we 
start in the deep-ferromagnetic region of the Heisenberg model. For the 
order parameter to jump from a finite value at gn = to a value of zero 
for any gn > there would have to be a first-order phase transition at the 
edge of the phase diagram. However, this edge is already well characterized, 
since it is the ferromagnetic phase of the Heisenberg model itself. Known 
properties of this phase are inconsistent with it being the location of a first- 
order phase transition in the same order parameter due to the following. 
A symmetry-breaking first-order phase transition is of the type associated 
with a tricritical point, which is controlled by a sixth-order Landau effective 
potential as shown in Fig. 2 (solid line) at the point of phase transition |14j . 
For gn > the two side-minima would hypothetically lift up, leaving only 
the minimum at zero order parameter. On the phase transition point, both 
the phase with zero order parameter and the two instances of the broken- 
symmetry phase exist in equilibrium. Such phase mixing would be easily 
observable through an un-sharp order parameter and a latent heat (range 
of internal energies). However, such a state simply does not exist in the 
Heisenberg model. In the deep ferromagnetic region its Landau effective 
potential is widely believed to look like the dashed line in Fig. 2 with only 
two minima which is incompatible with the tricritical behavior. 

Therefore, the hypothesis that the ferromagnetic phase does not enter 
the interior of the phase diagram seems to require behavior of the Heisenberg 
model inconsistent with known behavior. If a first-order phase transition is 
present it cannot be a conventional one. However, the argument above is 
somewhat heuristic and non-rigorous in that it involves Landau effective 
potentials. A more rigorous proof of the impossibility of a first-order phase 
transition as one enters the phase plane from the deep-ferromagnetic region 
of the Heisenberg model can be constructed as follows. Consider the corre- 
lation function C{q,r) =< ^tr(U^^4Ut^) >. In the ferromagnetic phase this 
shows spontaneous symmetry breaking through 

lim C{q,r) =< \m\ >^ (5) 

|g— r|— >oo 

where ^ 

^=yJ2^n (6) 

it 

where the sum is over individual 3-d hypersurfaces and V is the 3-volume. 
The expectation value above includes an average over different Euclidean 
times as well as gauge configurations. Although a finite lattice is used for 
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Figure 2: Landau effective potentiais at a symmetry breaking first-order 
transition (solid line) and in the deep-ferromagnetic region of a system with 
a higher-order transition (dotted line). Scales are arbitrary. 

defining quantities, the infinite lattice case is implicitly being considered 
here. 

Now let us calculate the derivative of this correlation function with re- 
spect to qh at qh = 0. Note that the calculation is performed in the 
Heisenberg model itself (earlier shown to be the gn = limit), but in order 
to determine what function to calculate the expectation value of, we must 
consider the gn 7^ theory. To do this, rewrite the horizontal U links in 
terms of gauge fields A, as in the usual calculation of the continuum limit: 



+l3vJ20---;;^'^i^^Pi^9HAfi^i ■r)Uji+i^4exp{-igHA^^i^ ••?)f^l,4)) (7) 



where i,j run from 1 to 3. No approximation has been made. For finite gn 
this is not a practical decomposition, since the A integrations still need to 
obey the Harr measure for the C/'s, however it is useful in the limit gn 0. 
Taking the limit gn — )• gives 




4+/3v.X(l-itr(Ug+j,4Utj) 



(8) 



n \i,j>i 
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where Fij is the abehan field strength tensor (the non-abehan part having 
a factor of gn)- This gives the multiple 3-d Heisenberg models as expected, 
but also a disconnected set of 3-d free field theories. The connected part of 



dC {q, r) 
dgH 



(9) 



is given by 



-/3y < tr(U5,4Ul Jtr((iAg,i • f)Ug+j,4Ut J > . 



n,i<4 



(10) 



and a similar term with the A in the other position. This is to be computed 
in the gn = action of Eq. [8] which is even in the transformation A — )■ 
—A. As a consequence, any expectation value containing a single factor 
of A is zero. Disconnected parts similarly vanish or cancel. This almost 
trivial observation means that the first derivative of the order parameter with 
respect to gy at = vanishes, and therefore exists. If the first derivative 
exists at a point, it follows from an elementary theorem of analysis that 
the function itself is continuous at that point [15]. However, if the order 
parameter is a continuous function of gy at = then it cannot drop 
suddenly to zero here in a first-order phase transition, which would require 
a discontinuity. Therefore, the ferromagnetic phase must enter the phase 
diagram. In fact, since gn does not affect the phase transition at lowest 
order, the line would be expected to enter at a 90° angle. 

Since, as argued before, the entire interior of the phase diagram is an 
SU(2) LGT if unequal lattice spacings are considered, this means that there 
must be a ferromagnetic phase in SU(2) as well. The SU(2) LGT has long- 
range order. The order parameter is the Coulomb-gauge magnetization. As 
argued before, due to the symmetry breaking nature of the phase transition, 
the phase transition line must continue to another edge, the only possibility 
being the lower right corner of Fig. 1. In doing so it clearly must cross the 
Pv = Ph line, which is the ordinary SU(2) LGT with equal lattice spacings. 
A phase transition in which the Coulomb Gauge remnant symmetry breaks 
spontaneously is known to be deconfining[3j. Therefore the zero physical 
temperature (infinite lattice) 4-d SU(2) lattice gauge theory must have a 
deconfining phase transition, contrary to the usual assumption. To prove 
for SU(3) or SU(N) one only has to replace the 0(4) Heisenberg model with 
the SU(N)xSU(N) spin model. Since these all have ferromagnetic phases, 
the argument goes through in the same way for them. 

This argument has similarities to another approach which had the same 
conclusion [16j. In the fundamental-adjoint plane of SU(2) lattice gauge the- 
ory with couplings /3f and /3a_, there is a well-known first-order phase transi- 
tion which starts at the 4-d Z2 LGT transition (at (3a = oo) and ends in the 
middle of the diagram at (/3i?, Pa) = (1-48, 0.90) [17j. This has been seen 
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as a critical point, below which an analytic path exists between the strong- 
coupling confining region with the weak-coupling continuum limit. However, 
that would require the first-order phase transition to be a non-symmetry- 
breaking transition. If it were symmetry-breaking, on the other hand, then 
the end of the first-order line would be a tricritical point, and the transition 
would have to continue as a higher-order one, bisecting the entire coupling 
plane into symmetry-broken and unbroken sectors, as above. Interestingly, 
a critical point and a tricritical point are rather easily distinguished through 
the scaling behavior of the attached first-order transition. In particular, the 
scaling relationship between the latent heat and the size of the metastability 
region is linear in the critical case and quadratic for the tricritical case[l6j. 
One simply monitors the shape of the growing hysteresis rectangle in the 
energy-coupling plane while moving up the first-order line. The numerical 
evidence points convincingly to the tricritical case[16]. This again implies a 
symmetry-breaking phase transition separates the strong and weak-coupling 
regions. One of the strengths of this energy-scaling argument is that it does 
not require identification of the symmetry or order parameter involved, and 
no gauge- fixing is required. The tricritical case simply requires some sym- 
metry to break. However, it seems likely it is the same symmetry breaking 
as studied above, something relatively easy to check. 

It is interesting to consider possible mechanisms which could drive 
the SU(2) transition. In the U(l) case, confinement arises from abelian 
monopoles, which are strong-coupling lattice artifacts [18j. The confined 
phase occurs when the monopoles form percolating loops. Confinement can 
be prevented by suppressing monopoles with a chemical potential [19]. Some 
time ago a gauge invariant monopole was introduced which could play the 
same role for the SU(2) theory[5j. It carries S0(3) and Z2 degrees of freedom 
which in some sense cancel each other, so was named an SO(3)-Z2 monopole. 
Suppressing these monopoles, together with a positive plaquette restriction, 
appears to prevent the transition to the confining phase [3]. The system stays 
in the spontaneously magnetized phase all the way to zero coupling when 
this constraint is imposed |20j. Lattices to 60^ have been measured and the 
lattice spacing is such that these should definitely be in the confining region 
if universality applies. It is interesting to monitor these monopoles in the 
standard Wilson theory on periodic lattices. They form percolating clusters 
in the crossover region where confinement occurs and de-percolate at weaker 
couplings. Fig. 3 shows the percentages of percolating lattices as a function 
of /3 for various lattice sizes. The finite-lattice percolation transition can 
be taken to occur at the 50% level. Also shown is the extrapolation of the 
percolation transition to the infinite lattice. The extrapolation shown uses 
the finite-lattice shift equation for a higher-order phase transition [21j 



using the expected value = 1.7 from Ref.[4J. Here L is the linear lattice 
size and /3cL is the apparent critical point on the finite lattice. A fit to the 



l/PcL = l//3c + cL 
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Figure 3: Percolation probabilities for SO(3)-Z2 monopoles as a function of 
/3 for various lattice sizes (diamond 16"^, * 18^, x 20^, box 2A^, triangle 30'', 
+ 40^). Also shown is extrapolation of 50% point to the infinite lattice (uses 
right y-axis). 



data for c and /3c gives /3c = 3.161(2). If instead, u is assumed to be un- 
known and determined from the fit then the best fit is obtained for v = 1.1 
yielding /3c = 3.193(1). Setting v = 0.5 gives /3c = 3.226(2). The uncertain- 
ties in /3c from the fixed-i' fits are much smaller than the differences in fits 
with different i^'s, so the majority of uncertainty in /3c is from the extrap- 
olation. Allowing a broad range for this exponent from 0.5 to 1.7 gives for 
the infinite lattice /3c = 3.19(3). It is very interesting that this is consistent 
with the position of the deconfining Coulomb-gauge magnetization transi- 
tion determined from data-collapse fits to scaling behavior, which yielded 
an infinite lattice critical point of 3.18(8) J^. It seems unlikely that these 
agree by mere coincidence. The percolation study was done after the other 
study was completed and released. Being that the percolation study was 
performed on the standard Wilson theory using conventional heat-bath and 
over-relaxation updates with periodic boundary conditions, the coincidence 
of these results for /3c gives additional confidence that the open boundary- 
condition Coulomb-gauge methods used in Ref. [4] are reliable. 

To summarize, it has long been known that lattice gauge theories have 
a lot in common with spin models. This becomes especially apparent in 
certain fixed gauges for which the remnant symmetry can be matched onto 
a spin model. The result presented above implies spin and gauge theories 
have even more in common, namely long range ferromagnetic order. This 
has important consequences for the non-abelian lattice gauge theories SU(2) 
and SU(3) which were previously assumed to be exceptions to this behavior. 
The relatively large critical exponent, ~ 1.7 for SU(2), explains how such 
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a transition could have been missed, since the corresponding specific heat 
singularity is very soft. 

The ordered continuum limit means that the quenched lattice gauge the- 
ory without fermions does not confine in the continuum limit, an idea which 
has been suggested previously[22J. Because good phenomenological evidence 
exists for quark confinement, however, some source of confinement must be 
sought, which need not necessarily be a linear potential. If light quarks 
are added to the theory, chiral symmetry is still expected to break sponta- 
neously through the formation of a quark condensate. This only requires 
a sufficiently strong force, not necessarily a confining one. If this collective 
state polarizes in such a way as to repel strong color fields, the chiral conden- 
sate could form a kind of bag surrounding mesons and baryons, contributing 
a confining-like term to the force over a limited range, though confinement 
would not be absolute [23]. Also, if a quark were to find itself a long distance 
from its partner antiquark, a quark/antiquark current in the chiral conden- 
sate could quickly generate local partners. This is somewhat different from 
generating quark anti-quark pairs from gluon "sparking" which originate at 
the same location, but that too can prevent isolated quarks from existing 
provided there is sufficient energy in the bond. Thermodynamics of the 
quark-gluon plasma is also modified in this picture, since confinement is no 
longer in the gluon sector. A phase transition or at least a rapid crossover 
to a chiral-symmetric phase at high physical temperatures will undoubtedly 
still exist, which, if chiral symmetry breaking is related to confinement as 
above, will also be deconfining. 
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